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Abstract 

We present an application of error theory using Dirichlet Forms in linear partial differential 
equations (LPDE). We study the transmission of an uncertainty on the terminal condition to 
the solution of the LPDE thanks to the decomposition of the solution on a wavelets basis. 

We analyze the basic properties and a particular class of LPDE where the wavelets bases 
show their powerful, the combination of error theory and wavelets basis justifies some hypothe- 
ses, helpful to simplify the computation. 

1 Introduction 

The study of the sensibility of the solution of a Partial Differential Equation to a small perturbation 
of the starting condition is a key problem in applied mathematics. The classical approach is to 
define a basis of "perturbation functions" and to evaluate the variation of the PDE's solution alongs 
each "direction"; this method is known as Gateaux derivative of the solution. But this approach 
assumes that the perturbation of the starting condition is deterministic, on the other hand the 
starting condition is often estimated by means of some measurements, i.e., with mathematical 
language, make a statistic; therefore the result of this estimation is a random variable, with a few 
known moments 1 . 

*I am grateful to Nicolas Bouleau for his helpful comments, work done partially during a visit to Institut Mittag- 
Lefiler (Djursholm, Sweden), 
email: simone.scotti@unito.it 

address: Via Real Collegio 30, 10024 Moncalieri (TO) ITALY 

1 Generally we know the mean, the variance and, maybe, the skewness and the kurtosis. 



The probabilistic nature of the uncertainty on each estimation pushes us to give up the Gateaux 
derivative approach; we follow the methodology suggested by Bouleau, see [3], this idea yields a 
representation of small perturbation coherent with the truncated expansion of the perturbated 
solution by a small random variable. If the estimation of starting conditions is good the related 
uncertainty is very small, we may neglect order higher than the second, that is we choose to 
work with Gaussian distributions. The error theory using Dirichlet forms has two operators, one 
describes the variance, called "carre du champ", other one the bias, called the generator of the 
semigroup. 

The analysis of the uncertainty transfer of the terminal function into the PDE's solution requires 
to specify a representation basis, we choose to work with a wavelet basis and we prove some useful 
properties in this case. 

Summarizing, we propose a new approach to study the impact of uncertainty on the solution 
of a Linear Partial Differential Equation due to an random imprecision on the starting condition; 
this method permits, first of all, the study of variance of the LPDE solution, i.e. his sensitivity; 
secondly we can estimate a covariance between the LPDE solution at two different points of time- 
space domain. 

The paper is organized as follows: In section 2, we present the error theory using Dirichlet Forms 
and we model, mathematically, the uncertainty on the starting condition. Section 3 is a survey of 
wavelets theory. In section 4 we study the solution of a LPDE using wavelet's decomposition and 
we present some particular cases where the wavelet's properties play a crucial role. In section 5 
we describe the evolution equation for the operators of error theory in the LPDE case. Section 6 
is dedicated to an introduction to possible applications in finance. The results are summarized in 
section 7. 

2 Error theory using Dirichlet forms 

To make the paper self-contained, we recall the principal ideas of error theory using Dirichlet forms, 
for reference see Bouleau [3]. 

The first studies of error propagation, owing to Legendre, Laplace and Gauss, date back at the 
beginning of the 19th century, show the probabilistic nature of errors. But the classical approach 
in the study of sensibility is to evaluate the derivative in the sense of Gateaux with respect to the 
erroneous parameters. Clearly this approach is coherent if the nature of incertitude is determinist 
(i.e. the "good" value exist) and it is incoherent with a probabilistic nature. 

The error theory using Dirichlet forms permits to study the sensibility with respect to a pa- 
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rameter, characterized by a random uncertainty. 

In order to introduce this theory, we analyze a particular case, we start with a parameter <tq 
and we want evaluate F(ao), where F is a nonlinear regular function. If the quantity ctq has a small 
centered error, we can model this parameter with a random variable a = <jq + eA[a] + e^[a]Y , 
where e is a small constant that defines the order of magnitude of the perturbation, A[a] represents 
the renormalized bias of the random variable with respect to the parameter value oo, j[<j] the 
variance and Y is a Gaussian standard random variable. 

If we study the bias and the variance of F(a), use of Taylor's formula show 

E[F{a) - F(<r )] = eF'(a )A[a} + ±eF"(<T )"f[a]+o{e) 
E[{F (a) - F(* )} 2 } = e[F'(a )] 2 -,[a}+o(e) 

We emphasize the following remark. 

Remark 1 The crucial difference between the deterministic approach and the probability one is in 
the bias term that depends on the variance of the random variable a if the function F is nonlinear; 
therefore this effect is said stochastic. 

The original idea of error theory using Dirichlet Forms is to stop the Taylor expansion at the 
first order in e and to forget the higher orders; and after search for a theory with two operators 
with the following chain rule: 

A[F[a)\ = F'{a)A[a\ + \F"{a)T[a] 
T[F(a)] = [F'(a)} 2 T[a] 

where T plays the role of a variance-covariance operator and A a bias operator. 

This theory exists and is the theory of semi-group, the operator A is the generator of the 
semi-group and the operator T is the carre du champ associated with the Dirichlet Form of the 
semi-group, see Albeverio [1] for a general introduction. 

We begin the introduction of error theory using Dirichlet Forms with the main definition. 

Definition 1 (Error Structure) An error structure is a term (CI, J 7 , P, D, T) where: 

1. (Cl, T , P) is a probability space; 

2. D is a dense sub-vector space of I? (fi, P); 

3. r is a positive symmetric bilinear application form BxD into L 1 (fi, J 7 , P), such that if F 
and G are two functions belonging C 1 and Lipschitzian class, u and v £ D, then we have 
F(u) and G(v) £ D and 
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T[F(u), G(v)} = F'(u)G'{v)T[u, v] Fa.s; 



4- The space D equipped with the norm \\u\\q = {\\u\\ L 2 +E[r[u, u}]) is complete. 

When the two arguments of gamma operator are the same, we can simplify the notation T[u] = 
T[u, u]. 

The knowledge of the carre du champ T and the probability P permits to associate a unique 
strongly continuous contraction semi-group (Pt) t > thanks to Hille Yosida theorem (see Albeverio 
pages 20-26). This semi-group has a generator (A) defined on a domain VA subspace of D and it 
is a self-adjoint operator that satisfies the chain rule (1) when F G C 2 , a G VA and 7(0-] G L 2 (P) 

The operator T is bilinear, like the variance-covariance one, that makes computations awkward 
to perform. It is possible to surmount this difficulty by means of a new operator, the sharp, which, 
in some sense, is a square root version of T, see Bouleau [3] page78-80. 



Definition 2 (Sharp operator) Let (fi, F, P, D, T) be and error structure and (^fi, F, PJ be a 
copy of the probability space. If the space D is separable, there exist an operator sharp, denoted 
()*, with these three properties: 



We illustrate the definition of error structure and related topics by means of an example. 

Example 1 (real space) The simplest structure on the real space is the Or stein- Uhlenbeck struc- 
ture: 



where B(M) is the Borel a -field ofR, fi is a gaussian measure and H 1 ^) is the first Sobolev space 
with respect to the measure fi, i.e. u G H 1 ^) if u G L 2 (fi) and v! in the distribution sense belongs 




1. Vm G D, u* G L 2 (P X P); 




3. Vu G B and F G C 1 , (F(u))* = F'{u)u*. 



(fi, F, P, D, T) = (M, B(R), h,H\h),T[u, u] = {u 1 } 2 ) 



to L 2 (ii). 



The associate generator has the following domain: 



T>A ={b£ L 2 (pL) : u" — x f' in the distribution sense belongs to L 2 (/x)} 



and the generator operator is 
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where I is the identity map on K 



The associate sharp operator is 



(u)*(x, y) 



u'(x)G(y) 



where G is a fixed function such that 



G 2 (y)d^y) = l 



R 



2.1 Error structure on functional spaces 

The idea of error structure also can be used when the probability space fl has infinite dimension, 
in this case some constraints must be verified see Bouleau [3] pages 59-64. In particular, afterwards 
in this paper, we need to define an error structure on the space of functions L 2 . 

Let fl be a space of functions, we assume that it exists a basis (perhaps numerable) {V'ijieN of 
this space; therefore if / S then 



The simplest way in order to define an error structure on the space f2 is to randomize the 
coefficients a,, to define an error structure on each subspace generated by each element of the 
basis and take the product structure as the error structure on the functional space fi. Bouleau 
has proved that an infinite product of error structures is still one, under some constraints on the 
r operator, see [3] page 62. 

An easy way to force the wellness of the product of error structures is to truncate the previous 
sum, see Bouleau [3] page 84 for a practical example. 

On each factor we have the choice on the structure of gamma operator; a first idea is to suppose 
that the uncertainties on two different factors are uncorrelated: 

Hypothesis 1 (independence) The operator T verifies: 



This hypothesis simplifies the course of the analysis of uncertainty diffusion, we emphasize the 
nature of this hypothesis in the following remark. 
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F[a m , a n ] = when m ^ n 



(2) 



Remark 2 It is clear, first, that the hypothesis of uncorrelation must be have an external jus- 
tification, the error theory using Dirichlet forms cannot force a particular structure on the error 
correlation. Secondly, the choice of the basis is crucial, since, if we start from an error structure, 
associated at a basis, that verifies the hypothesis (1), and we change the basis, the image error 
structure through the new basis presents some correlations except in exception circumstances. 

Therefore, it is clever to choice a basis with some features in order to use these properties to 
justify the hypothesis of uncorrelation. 

Lastly we have to specify the choice about the operator T on each subspace generated by each 
element of the basis; many options are possible but we prefer the following 

Hypothesis 2 (proportionality) The operator T verifies: 

T[a n , a n ] = 7(71) a 2 (3) 

The advantage of this hypothesis is to yield a constant proportional error, i.e. each factor on 
the basis decomposition contributes towards the uncertainty proportionally to their contribution 
at the decomposition of the function; therefore, if a function as not projection on a subspace, the 
same subspace do not contribute to the uncertainty. The function 7(71) permits to leave some 
freedom at the class of error structure. 

Considering the fact that we have an independent error structure on each subspace generated 
by the basis decomposition, thus the Mokobodzki hypothesis is verified, the error structure admits 
a sharp operator, see Bouleau [3] pages 78-82, and it verifies the properties: 



<4 = \/j{n)a n a n 

( 4 ) 

Ep£] = 1 



3 Wavelets theory 



In this section we give a short presentation of the theory of wavelets, in order to preserve the 
self-contained of this paper, this introduction follows Hardle et al. [10] and Mallat [12]. We fix 
a function, called father wavelet, G L 2 (R), such that the family {4>o,k = <j>{- —k),k G Z} is an 
orthonormal system 2 . 

We define the linear space (sub-space of L 2 ) 



Vo = /(•) 



/(•) = c ^,k{-) 



The easy way is to consider a compactly supported 



From this original space we can define a chain of sub-spaces {Vi}i£% by the relation 



/(•) e V iff /(2.) g K-i 



these spaces are called "generated" by the function (j>. Mallat and Meyer introduce the concept 
of multiresolution analysis in the years 1988-1990 

Definition 3 (MRA) A chain of sub-spaces {Vi}i^z, "generated" by a function (f> is called a 
Multiresolution Analysis [MRA] if the two following conditions are held: 

• Vi C Vi+i for all i and 

• \J i Vi is dense in L 2 

In this case the function <f> is called the "father Wavelet". 

To define an orthogonal basis of the L 2 space we must define a sequence of orthogonal spaces, 
since the chain {Vi} is decreasing sequence. Define Wi the orthogonal complement of Vi into V$_i, 
for all i; we find that Vi = Vo ffi ©}=i Wj and, thanks to the second property of MRA, we have: 



We can fix an orthonormal basis {ipik}keZ m eacn space Wi, Mallat and Meyer show that we 
can fix ipi t k{-) = v2~^V , (2 - '? • —k) where tp(-) is a function depending of (/>(•), see Mallat [12], page 
233 for the explicit relation; the function ip(-) is called the "mother Wavelet". 

We conclude that the original function f(x) has a unique representation in term of the following 
series: 



3.1 Daubechies wavelets 

The construction of the wavelet basis depends on the choice of the father wavelet; this choice is 
not constrained, many options are possible, in this section we present a particular class of wavelets 
called Daubechies Wavelets, see Daubechies [5]. 

The first historical wavelet is the Haar's basis, its father wavelet is (f>(x) = l[o, i]l this basis 
has some advantage, in particular the boundary support and the quick computation; unluckily the 
Haar mother wavelet is discontinuous, therefore this basis misapproximates all continuous function, 
this roughly approximation induces that the coefficients k do not decrease fast with the rescaling 
index i. 



oo 



3 = 1 




k 



k 
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Another possibility is the Riesz's class, but the choice of the Riesz's bases approach cause 
uncompactly supported father and mother wavelets, see [10] chapter 6, and it is clear that a 
function with compact support is easy to treated numerically. 

Actually the simpler known class of wavelets with continuous and compactly supported father 
and mother wavelets is the Daubechies' wavelets, see [10] chapter 7 or [12] pages 246-251. 

The definition of Daubechies' father wavelet starts from its Laplace transform. 

oc 

0(w) = JJmo(2- J 'u>) 

3=0 

where m is a 27r-periodic function, that must verify some constraints, see [10] chapter 5 and 6 or 
[12] pages 222-231; in particular an usual choice for m is a trigonometric polynomial: 

k 

The Daubechies' basis of order p must verify that the function m (cu) has a zero of order p for 
u = tt; we conclude with the equation of the Fourier transform of mother wavelet: 

$(u) = -±=m (2-W7r) e -*t0 (|) 
the principal advantages of these bases are: 

• regularity: the Daubechies' father and mother wavelets are uniformally 7-Lipchitz where the 
parameter 7 grow with the order p, in particular if p is grater than three the wavelets are 
differentiable, see Daubechies et al. [7]. 

• compact support: the support of Daubechies' father and mother wavelets are length 2p — 1. 

In conclusion we have a representation of the function f(x) via a Daubechies' decomposition 
on wavelets: 



/(*) = « fe < k (*)+EE Ik k ip° k (x) (5) 

k i k 

4 Linear partial differential equation 

In this section we want study the interaction between wavelets and linear partial differential equa- 
tions with a terminal condition f(x). In the previous paragraph we have decomposed the function 
f(x) into a Daubechies' wavelets basis; now we study the evolution of this decomposition through 
a linear partial differential equation and we focus on the properties of this evolution. 
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We consider a Linear Partial Differential Equation 




a 2 (x, t) d 2 Q 
2 ~dx 2 



{t,x)+n(x,t)^- = 



< 



(6) 



Q(T, x) 



The Feynman-Kac formula tells us that the solution of LPDE (6) is 



Q(i, x)=E[f(X T )\X t =x] 



(7) 



where the Ito process X t verifies the following SDE 



dX t = n(t, X t )dt + a{t 7 X t )dW t 



(8) 



With Wt is a Brownian motion. 

Since the linearity of the PDE (6), the terminal condition and the solution admits a decom- 
position on basis, in particular a wavelets basis. Then, if the terminal condition is written as 
following 



we can study the evolution of each factor and prove that the solution has a similar decomposition. 

Notation 1 In order to simplify the notation we shift the time coordinate in such a way that we 
bring the final time T at zero. Clearly each time t smaller than T became negative and the terminal 
condition became Q(0, x) = f{x). 

4.1 Diffusionlets 

In this subsection we analyze a particular case in which the wavelets properties have a key role; 
a basilar study has been done by Shen and Strang, see [15], when the PDE is the heat equation 

U t (t, x) = cAU(t, x). 

Shen and Strang introduce the notation of mother and father heatlets (given a choice of wavelets 
basis), these are the heat evolution of the mother and father wavelets. 

In a similar way, if we consider a LPDE and a wavelets basis, we can define an object called 
Diffusionlet: 



/(*) = + E E a, ki% k (?) 



k i k 



Definition 4 (Diffusionlets) The solution of a LPDE with a wavelet as terminal condition is 
called diffusionlet. 
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The diffusionlets associated at a LPDE are denoted: 



we can make some remarks: 



cf>° k (X )\X t =x 
ipP k (X )\X t =x 



(9) 



Remark 3 In Daubechies' case, the father and mother wavelets are compactly supported and 
bounded, therefore the variance of the solutions Q D and fy D are, generally, smaller than the vari- 
ance of the solution Q(t, x) 

Remark 4 The Daubechies' mother wavelet %l) D {x) and its rescaling functions ipf k (x) have p 
vanishing moments (where p is the order of Daubechies 's wavelet), so the associate diffusions depend 
mainly on the assymmetries and we can suppose that the contributions of high order wavelets vanish 
very quickly with the time. 

Remark 5 The study of the solution of the LPDE (9) requires to solve the same LPDE with each 
wavelet basis as terminal value. 

The least remark underline the main problem with this approach, we have earned a high pre- 
cision on the estimation of the solution but the price to pay is that we need to solve many times 
the same problem in order to have a good estimation. But we introduce a class of LPDE where 
this difficulty has an easy answer. 

We consider a special case of the previous LPDE, we fix the function fi(t, x) = rx and the 
function a(t, x) to be a power of the variable x. 

The LPDE becomes 



f dQ a 2 ^ d*Q dQ 



(10) 



Q(0,x) = f(x) 



and the associate diffusion is 



dX t = rX t dt + aX^dWt (11) 

A change of numeraire proves that the good solution of LPDE (10) can be write as Q(t, x) = 
Q(t, x)e~ rt where Q(t, x) verifies the following LPDE 
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f dQ a 2 x 2X d 2 Q 

(12) 

Q(0,x) = /(*) 

Remark 6 (Self-similarity) TTiis equation presents a reseating invariance: if Q(t, x) is a solu- 
tion of LPDE (12) with terminal value g{x), then also Q(a 2 ~ 2X t, ax) is a solution of LPDE with 
terminal value g(ax). 

We can prove some properties and theorems of the wavelets diffusion through a LPDE of type 
(10), each results is a generalization of an equivalent theorem on heatlets decomposition, see Shen 
and Strang [15], the crucial hypothesis used in all theorems is the linearity of PDE. 

Proposition 1 (Refinement) Let 4>{x) and il)(x) be, respectively, the father and mother wavelets 
of a wavelets basis. Suppose that the LPDE (12) has a unique solution for each terminal func- 
tion f(x); let $> D (t, x) and v E' D (t, x) be the solutions of the LPDE (12) with terminal conditions, 
respectively, the father and the mother wavelets. If 

4>(x) — h n (f)(2x — n) 

n G Z ( 13 j 

ip(x) = 2 ^2 9n H 2x - n ) 

with 

nez 

Then 

$ D (t,x) = 2 ]T h n <S> D ((2( 2 - 2A ^, 2x - n) 

nez 

(14) 

V D (t, x) = 2^ 9 „$ D (2( 2 - M 't,2 I -r i ) 

n6Z 

Proof. It is easy to check that <& D (2 {2 ~ 2X) t, 2x - n) is the solution of LPDE (12) with terminal 
value 4>(2x — n), see remark 6. Then the evolution of the functions <j)(x) and ?/>(x) given the left- 
hand sides of equations (14) and the decomposition of the same functions, thanks to relations (13) 
given the right-hand sides; by uniqueness the two solutions must be equal. 

□ 

Theorem 1 (Diffusionlets decomposition) Suppose that f(x) belongs to L 2 and then f{x) ad- 
mits a decomposition of type (5) and suppose that LPDE (12) admits a unique solution when the 
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terminal condition is Q(0, x) = f(x) then the solution of LPDE (12) with terminal value f{x) is 
given by 

Q(t, x) =J2ak<S> D (t,x-k) + J2J2 k * D (2 (2 " 2A)I t, V x-k) (15) 

k i k 

Proof: Since f(x) g L 2 (R) and |</)jp V^fcW} is an orthonormal basis of L 2 (R), the wavelets 
expansion of f(x) converges to f(x) in L 2 -norm. 

It is easy to check that * D (V 2 ~ 2A)i t, V x - kj is the evolution of ^ k (x) and $ D (i, x-k) is 
the evolution of ^ k {x), see remark 6. 

By uniqueness of the solution the diffusion of the wavelets expansion of terminal condition f(x) 
converges to the diffusion of f(x) in L 2 -norm 

□ 

Remark 7 (Diffusionlets Advantages) The key advantage of the diffusionlets is the indepen- 
dence of the initial state. Therefore we can store the solution of the LPDE (12) with the father 
wavelet as terminal condition, this solution can be estimate with an high degree of precision due to 
the compactly supported and bounded father wavelet in Daubechies case. 

The solution of the LPDE (12) with terminal value f{x) can be evaluate using this strategy: 

1. study the decomposition coefficients of the function f(x) into the wavelets basis, using FWT 
(fast wavelet transform) as an example, see Mallat [13]; 

2. reconstruct the solution with f{x) as terminal condition using the store solution for the father 
wavelet, the coefficients estimated and the result of theorem 1. 

The second, crucial, advantage is the fact that Daubechies wavelets of order p have p vanishing 
moments, this fact, combined with the smoothing property of parabolic partial differential equation 
(property true for LPDE (12) far from x = 0) force a fast convergence of the mother wavelet to 
with time t; this effect is magnified by the scaling effect when A < 1, see equation (15). 

Remark 8 (Diffusionlets disadvantage) Due to the linearity of the LPDE (12), diffusionlets 
are not compactly supported when t ^ 0. This fact force the necessity to estimate an "essential 
support " for the diffusion of father and mother wavelets, i. e. the region where the norm of father 
and mother diffusionlets are bigger than a reference level e, related to the asked sensibility of the 
searched solution. 

This "essential support" given a length, depending on time t, useful to determine the number 
of decomposition elements for each scaling level, i.e. the number of k that it is necessary consider 
for a well-estimation of the solution before changing the coefficient i. 
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The smoothing effect, emphasized between the advantages, permits to define an order /(e) 
beyond which the remainder is smaller than e 
The approximated solution becomes: 

K (e) 1(e) Ki{e) 

Q{t,x)~ Uk$ D (t,x-k)+J2 E A,fc^(2 (2 - 2A H ?x~k) (16) 

fc=-A' (e) i=0 h=-Ki(e) 

We conclude with a negative remarque on the diffusionlets basis: 

Remark 9 (non- orthogonality of diffusionlets basis) The class of functions generated by the 
diffusion of a wavelets basis is a basis of the L 2 -subspace characterized by the diffusion itself , i.e., 
for each time s, the smaller subspace of L 2 that contains the functions that can be a solution of the 
LPDE at time s with a terminal value belongs to L 2 . 

However, the basis {$f (t, x), $£ k (t, x)} . , where <&f (t, x) = $ D (t, x - k) and ^ k {t, x) = 
ty D ^2^ 2 ~ 2A ) 1 1, 2 l x — kj , is not, generally, orthogonal, especially owing to the scaling factor in 
time. 

5 Sensibility of LPDE solution 

In this section we suppose that the terminal condition f{x) of a LPDE are erroneous and we study 
the diffusion of this uncertainty. The starting point is to define an error structure on a functional 
space, see subsection 2.1. We use the decomposition of f(x) into a wavelets basis, see equation (5); 
and we set the coefficients au and k to be random. 

We define an error structure on each subspace, generated by each element of the wavelets basis, 
in accord with hypotheses 1 and 2, we can assume that the error structures on each subspace 
are independent and the uncertainty is proportional to the estimate parameter, i.e. au or k 
depending on the cases. Now we can study the variance caused by the uncertainty on the terminal 
value. 

5.1 uncertainty on the solution 

We can prove three results: 

Proposition 2 (Variance of terminal condition) The terminal condition f(x) has the follow- 
ing variance: 

3 The proof of this fact is very simple, we must check only the independence of each vector of the basis, that 
follows on the uniqueness of the solution. 
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r [/(*)] = £>(fc)<4 [V&M] 2 (17) 



I 2 

k i k 

Proof. We start with the computation of the sharp of the function f(x), the definition of the 
sharp operator (2) and the chosen structure (4) give the following relation: 



k i k 

(18) 

= v'tW & k &k<Po,k( x ) + Vt(M0 P%,kPi,k i>f,u ( x ) 

k i k 

Now the second property of operator sharp (see definition 2) gives the result (17). 
The previous theorem as an equivalent for the solution of LPDE. 

□ 

Theorem 2 (Variance of LPDE solution) The solution of LPDE has the following variance: 

r[Q(t,x)] =$>(*W [< fc (^)]%£E^ fc )#U Mk^ x )\ 2 ( 19 ) 

k i k 

The proof is similar to the previous theorem, we start with the sharp decomposition: 



k z k 

= J2 ^W) otk s k x) + E fc ) A, fc ft, k 



(20) 



i k 



Now the second property of operator sharp, see definition (2), gives the result (19). 



□ 



Remark 10 This uncertainty is easy to estimate due to the fact that of diffusionlets are universal, 
in the sense that they are independent of the initial state, see remark 7. 

Theorem 3 (Covariance of LPDE solution) The solution of LPDE has the following covari- 
ance: 
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r Q(t, x), Q(s, y) 



J2^ a l x) $£ fc (s, y) 



k 



(21) 



i ft 



The proof is equal to the previous theorem. 



□ 



6 Applications in Finance 

In the domain of mathematical finance, the Partial Differential Equations have a key role. Excepted 
few cases, the equations in a financial model do not have a closed form solution, therefore a 
numerical approach is mandatory. The Feynman-Kac formula gives us a bridge between the PDE 
and SDE, hence it is possible to choose the numerical method, PDE or Monte-Carlo, depending 
on the time-efficiency; generally PDE approach is used in the case of a low number of variables. 
A first application is the Cox-Ingersoll-Ross stochastic differential equation, see [4]: 



This case has a relative interest, since an exact characterization of the solution exists, see 
Lamberton and Lapeyre [11] or Shreve [16]. The principal advantage consists on the possibility to 
test the wavelets procedure. 

A second, more relevant, application is a simplified SABR model; Hagan et al, see [9], emphasize 
the incoherence of the dynamic behavior of log-normal model, the well-known Black Scholes model, 
compared to the behavior observed in the marketplace. In order to eliminate this problem, Hagan 
an Woodward, see [8], propose a local volatility model in which the forward value satisfies 



where A is a fixed parameter, that takes values between and 1, estimated on the market. This 




when the mean a is zero; the associate PDE is 




d 2 U 



dx 2 



dF t 



at F t x dW t 
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model is the starting point for the SABR model, the SDE is of type (8); therefore the procedure 
described in this article can optimize the procedure of option pricing. 

Hagan et al, see [9], emphasize an other relevant aspect, actually market smiles are managed 
using Dupire's local volatility models, but a local volatility function different to a power 4 introduce 
an intrinsic "length scale" for the forward price, this inhomogeneous has an hard financial expla- 
nation. Therefore the model proposed by Hagan and Woodward has the principal problems of all 
local volatility model, i.e. a poor equivalent volatility for basket options and an unsmiled forward 
volatility, we propose to solve these inconveniences using the same strategy proposed by Scotti, 
see [14], that is consider a perturbation on the volatility parameter at and handle this uncertainty 
through an error structure; this study will be investigate in a following paper. 

7 Conclusion 

In this paper we have investigate the relations between three objects, i.e. the linear partial differ- 
ential equations, the error theory using Dirichlet forms and the wavelets. 

These three objects have a different origin and an application fields separate to date; this article 
shows as capitalize the advantages of wavelets bases in order to solve a LPDE and as exploit the 
wavelets as a decomposition basis to study the sensibility of the LPDE. 

In a particular case, when the LPDE has the following form: 

U t = a 2 x 2X U xx 

We have proof that the properties of wavelets are partly preserved, especially the invariance 
under a scaling and a translation; these properties permit a fast processing of the general solution 
of this type of PDE. The consider LPDE is the new key element of the non-lognormal financial 
models, studded by Hagan et al, see [8] and [9]. 

The principal new feature is the study of the sensibility of the PDE solution using error theory 
by means Dirichlet forms, introduced by Bouleau, see [3], this methodology has our permitted 
an evaluation of the sensibility of the solution w.r.t. an uncertainty on the terminal value, this 
inaccuracy on the payoff can model an imprecision on the final spot value, due to the spread bid-ask 
as an example. 

This paper is just a starting point for the study of inaccuracy on terminal value and non- 
lognormal models, in following papers will be evaluated the impacts of an uncertainty on the 
volatility parameter a and on the implied curve of volatility. 

4 In the case of a power the local volatility model is reduced at a simplified SABR model. 
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